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Abstract 



We study the decay of charged scalar and spinor fields around Reissner-Nordstrom 
(N ■ 

J> ■ black holes in de Sitter spacetime through calculations of quasinormal frequencies of the 

rn ; 

I fields. The influence of the parameters of the black hole (charge, mass), of the decaying 

^ . fields (charge, spin), and of the spacetime (cosmological constant) on the decay is analyzed. 

' The analytic formula for calculation quasinormal frequencies for a large multipole 



number (eikonal approximation) is derived both for the spinor and scalar fields. 

1 Introduction 

It is well known that dynamical evolution of field perturbation on a black hole background can 
be conventionally divided into three stages The first stage is an initial wave burst coming 
directly from the source of perturbation and is dependent on the initial form of the original field 
perturbation. The characteristic feature of the second stage is the damped oscillations whose 
frequencies and damping times are defined by the structure of the background spacetime and 
are independent of the initial perturbation. This stage can be accurately described in terms 
of the discrete set of quasinormal modes (QNM). And the third stage is an asymptotic tail 
behavior of the waves at very late time which is caused by backscattering of the gravitational 
field. 

The original interest to the study of QNMs is arising from the possibility to observe quasinor- 
mal ringing with the use of gravitational wave detectors as it follows from theoretical predictions 
(for review see, e.g., and references therein). Recently the interest to the study of QNMs 
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has been reinforced in connection with their possible relation to the thermodynamic properties 
of black holes in loop quantum gravity 

A new application of the quasinormal mode spectrum has also arisen from superstring theory 
iniinilZ!- According to the AdS/CFT correspondence, a large static black hole in asymptotically 
AdS spacetime corresponds to a thermal state in CFT, and the decay of the test field in the 
black hole spacetime corresponds to the decay of the perturbed state in CFT, such that the 
quasinormal frequencies define the thermalization time scale [8]. Therefore, many authors 
focused their attention on the studies of QNMs for different asymptotically AdS black holes 

iHiiniiiniiiiiiiiiisi- 

There exist observational evidences that the universe is described by the general relativity 
equations with the positive cosmological constant ^1 |^ ^j. That observation attracted 
considerable researches' attention to the study of QNMs in asymptotically de Sitter spacetimes 

ii7i[iHiiinii2niinii22ii2ii2ai23. 

In ^3 the calculation of the quasinormal frequencies for the gravitational perturbations 
of the Schwarzschild de Sitter (SdS) black hole was carried out. In ^22jj the lower overtones 
of QNMs for higher dimensional SdS black holes were calculated. In the total spectrum 
of QNMs was obtained for the SdS black holes by numerical calculations and the excellent 
coincidence with the 6th order WKB method was shown. In [23^ the low-laying quasinormal 
frequencies of the SdS black hole for fields of different spin were calculated by the use the 6th or- 
der WKB and the Poshl- Teller potential approximations. In ^H] the authors derived an explicit 
expression for calculation of QNMs for the case of a near extremal SdS black hole. In [20] that 
expression was generalized to near extremal higher dimensional SdS and Reissner-Nordstrom 
de Sitter (RNdS) black holes. In [21] an analytical method was developed to study the quasi- 
normal mode spectrum of SdS black holes for the scalar, electromagnetic, and gravitational 
fields in the limit of nearly equal black hole and cosmological radii. 

The QNMs of the massless uncharged Dirac fields for the RNdS black hole are studied in 
PH] using the Poshl- Teller potential approximation. It was found that the magnitude of the 
imaginary part of the quasinormal frequencies decreases as the cosmological constant or the 
orbital angular momentum increases, but it increases as the charge or the overtone number 
increases. We note although that first the Dirac quasinormal modes were evaluated in |2S1 for 
Schwarzschild black hole spacetimes. 

Perturbations of a charged massless and massive scalar field was by calculation of its QNMs 
in the Reissner-Nordstrom (RN), were investigated in |l22j and jlHj, respectively. It was found 
that the neutral perturbations dominate at the stage of the "final ringdown" . In IHO] the 
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spin 1/2 Dirac particles with the positive, negative, and zero charge e in the presence of the RN 
black holes with charge Q was investigated. It was demonstrated that in late times the neutral 
perturbations dominate when eQ > and the charged perturbations dominate when eQ < 0. 

The decay of fields, which interact electromagnetically with the charge of a black hole in de 
Sitter spacetime, was not considered till now. The objective of the present paper is to study 
the decay of charged fields of different spin in RNdS background and analyze how the decay 
is influenced by a total set of parameters: the charge Q of the black hole, the cosmological 
constant A, and the charge e of the decaying field. 

The paper is organized as follows. In Sect. 2 we consider the spinor field decay, in Sect. 3 
we consider the scalar field decay, and in Sect. 4 we summarize the results obtained. 

2 Decay of the spinor field 

In Schwarzschild coordinates, the metric for the Reissner-Nordstrom de Sitter black hole can 
be expressed as follows 

ds^ = -fde + f-^dr^ + r^ {de^ + sin^ Odif^) (1) 



2M g2 ^ 



-2 3 



(2) 



where the parameters M, Q are the mass and the charge of the black hole, respectively, and 
A is the cosmological constant. We introduce the tortoise coordinate = J f^^dr. Then we 
write the metric function (jSJ in the form 

f = :r^{r-r^){r- r+) (r^ - r) (r - n) , 
3r^ 

where is the cosmological horizon, r_ is the inner event horizon, and r+ is the outer event 
horizon {vb can be found through r_, r+, Tc). One can find that the tortoise coordinate r* can 
be expressed as follows 



— In (r — r_) H In (r — r_| 



(3) 



In (rc — r) H In (r — rj,) 



Idl 
2dr 



, a = (-,+,c,6) . 

'a 

From (0) one can see that —>■ oo as r —>■ and r^, — > — oo as r 
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The wave equation of the massless Dirac field can be written as 



{d, + V, + eA,)] = (4) 

where 7" are the Dirac matrices, is the inverse of the tetrad , = | [7*^,7^] ^a^bv;fi is 
the spin connection. We take the tetrad as [SI] 

e^ = diag{f'\r"\r,rsme) . (5) 

We define $ = then using the ansatz 



where 



~ If ■ - / -u 1 

-jm I ^ /7^y^m+l/2 ) tOr J - i + ^ 



^i±i-' for 7 = / - - 

^ ^' 2i+2 '^l I 



we find for both signs of F and G, 



drl 

where 



+ (cu^ -V2)G = 



,,, = ±^ + ».^W = M^-^^2, (9) 

ar^ r r 

k=j + ^, j = / + ^ for Vi , 

A; = - (^J + , J = ^ - ^ for V^2 • 

The quasinormal modes are defined as solutions of (|7j), (jH)) satisfying the boundary conditions 

F in) ~ A±e±'"'-% G in) ~ 5±e±'"'-% ^ ±00 , (10) 

supposing Recj > 0, that corresponds to purely in- going waves at the event horizon and purely 
out-going waves the cosmological horizon. 

We notice that the potentials Vi, V2 give the same quasinormal frequencies jSH], [HO], there- 
fore in what follows we focus the attention on Vi. The effective potential Vi is a smooth 
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Figure 1: QN frequencies of the spinor field; red: e = 0.1, blue: e = —0.1, green: e = 0. 

-Imo) 




Figure 2: Imaginary part of on of the spinor field; red: e = 0.1, blue: e = —0.1, green: e = 0. 

-Imo) 




Figure 3: QNM of the spinor field parametrized by Q for A = 0, 0.2, 0.4, 0.6, 0.8, 0.1; red: 
e = 0.1, blue: e = —0.1, green: e = 0. 

function of r^,. It tends to constant values at the cosmological horizon and at the outer event 
horizon and has a maximum near the outer event horizon that allows to use the WKB method 
to calculate quasinormal modes. To compute QNMs from Eqs. ((Tj), we apply the WKB 
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Figure 4: QNM of the spinor field parametrized by A for e = 0.1, Q = 0, 0.2, 0.4, 0.6, 0.8. 



Figure 5: QNM of the spinor field parametrized by A for e = —0.1, Q = 0, 0.2, 0.4, 0.6, 0.8. 

method of the sixth order jH2j . 

We present graphically the results for the fixed overtone number n = and = 1 at Figs. 
1-5 for the charge of the spinor field e = 0, ±0.1. It is shown the dependence of QN frequencies 
(Fig. 1) and their damping rates (Fig. 2) on the charge of the black hole Q for different 
values of the cosmological constant A. On Figs. 1, 2 the dependences are shown for three 
values of the charge of the spinor field e = 0, ±0.1, which on the figures looks like a triplet of 
curves beginning at the same point at Q = 0. The curves of different densities on the graphs 
correspond to different e. On Fig. 1 the upper curve in the triplets corresponds to e = —0.1, 
the lower curve in the triplets corresponds to e = 0.1, the middle curve corresponds to e = 0. 
On Fig. 2 the lower curve in the triplets corresponds to e = —0.1, the upper curve in the 
triplets corresponds to e = 0.1, the middle curve corresponds to e = 0. The QN frequencies 
are monotonically growing with Q for e = 0, —0.1, although for e = 0.1 the behavior of the 
curves is qualitatively different. The curves for e = 0.1 have a minimum at some point Qq. 



-Imo) 
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The value of Qq is shifting to less values as A grows. Besides, Fig. 1 demonstrates that the 
QN frequencies become lower as A increasing. The position of a curve in the triplet says that 
the QN frequencies is lower for the field charge e > and higher for e < (Q > 0). The 
curves displaying QNM damping rate are similar for all A and e. They have a maximum at 
Qmax ~ 0.9. Until reaching Qmax the curves are slowly growing and after reaching Qmax they 
are sharply falling down. The greater A is the greater the rise of the curves is until reaching 
<3max- The value of Qmax is shifting to greater values as A increasing. The position of a curve 
in the triplet says that the decay of the field is slower for the field charge e < and faster for 
e > (Q > 0). The decay of the spinor fields becomes slower as A grows. Since the charges e, Q 
appear in Vi only as e^, or eQ one can also conclude that for Q < the conclusions become 
upside-down. These observations generalize results the work [30j obtained for asymptotically 
flat spacetimes. The peculiarities of the dependence of Imu, Keu on Q are more obvious when 
depicted on diagrams of the dependence of Imu; on Rea; and parametrized by the black hole 
charge Q. These curious pictures are presented on Fig. 3. 

The dependence of Imu) on Reu; is parametrized by the cosmological constant A, which is 
changing along the curves, is depicted on Figs. 4, 5 for the fixed black hole charge Q. The 
curves are approximately straight lines. This observation generalizes conclusions the research 
|25] obtained for the uncharged spinor field. Both the QN frequencies and their damping rates 
are decreasing with A, although the manner of the decrease for Recu and Imu is different. 

Making use of the first order WKB method we also found the asymptotic formula of large 
k for calculation quasinormal modes uj, 



where 



uj = Cok-iin + 1/2) Cf + C'/ + O (l/k) (11) 



Cf = -i- [-42Q^ + 3QVo (42M - 15ro + 2Arl) + 
rl (QOM^ - 60Mro + 9rl + QMAr^ - Ar^)] , 

yl/2 



4rf roCo'^' + (ro - 3M) 



ACo'^' {Q' + eQrl + rf (3M - 3ro + 2rlA)) 
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^1 — 2 ° ° 

(-30g^ + 90MgVo - 63MV^ - 33gV^ 
-18egV^ + 42Mr^ + SOeMQr^ - 6r^ 
+12egr^ + egV^A - 6Mr^A + r|^A + 2egr^A) 
(-42g^ + 126MgVo - 90MV^ - 45gV2 
+60Mr|j - 9r^ + 6gV^A - 6Mr^A + r^A) . 

By setting Q = one can make sure that the formula (fTT|) reproduces the corresponding formula 
of the paper obtained for the SdS black hole. We note that Imcu does not depend on the 
charge of the spinor field e in the limit of large k. Therefore the damping rate is equal for the 
arbitrary spinor field charge for large enough multipole moment. The dependence of Keu on 
k is linear for large k. The slope of the line on the plane (Reu;, k) is defined by the coefficient 
Co- These conclusions are confirmed by the results of [30] where the particular case A = was 
considered. 

Re CO 




Figure 6: QN frequencies of the scalar field; red: e = 0.1, blue: e = 0.05. 



3 Decay of the scalar field 

The wave equation of the complex scalar field has the form 

(P.ab - leAag^' {2<p.b - teAi,<P) - leAa.bg'"'' = (14) 

Here we choose the electromagnetic potential as At = —Q/r. Then, decomposing the wave 
function into spherical harmonics = ^im'^''m(t,r)Y'-m{0,ip) /r, the wave equation for 
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Figure 7: Imaginary part of uo of the scalar field; red: e = 0.1, blue: e = 0.05. 
-lma> 




Figure 8: QNM of the scalar field parametrized by A for Q = 0, 0.2, 0.4, 0.6, 0.8; e = 0.1. 
each multipole moment takes the form 

d'^u + 2ie—dtu - d^u + Vu = (15) 



\ r r'^ 6 J 

where is the tortoise coordinate and u {t, r) = 6"*"^*^ (r). Therefore we come to the equation 
for u (r), 

a?.M+ [i^'-K//]« = o, (17) 

Supposing Reu > solutions of (fTTj) have to satisfy the boundary conditions 

« (r,) ~ C±e±'"'-% r, ^ ±oo . (18) 
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To calculate quasinormal modes we apply the WKB method of the third order in the 
scalar field case. The motivation about the use of WKB method is the same as that in Sect. 2. 
The use of the calculation method of higher order |22i is overly cumbersome as in the case of 
scalar field Ve// includes uj so that it is necessary to solve simultaneously the equation for the 
search of an extremum and the equations of the WKB method itself. 

We present graphically the results for the fixed overtone number n = and multipole 
number / = 1 on Figs. 6-8. It is shown the dependence of QN frequencies (Fig. 6) and their 
damping rates (Fig. 7) on the charge of the black hole Q for different values of the cosmological 
constant A. On Figs. 6, 7 the dependences are shown for two values of the charge of the scalar 
field e = 0.1, 0.05, which on the figures looks like a couple of curves beginning at the same 
point at Q = 0. The upper curve in the couples corresponds to e = 0.1, the lower curve lower 
one of the couples corresponds to e = 0.05. 

The QN frequencies monotonically grow with Q and become lower as A increasing. The 
QN frequencies are higher for the fields with greater charge e. The curves displaying the QNM 
damping rates have a pronounced maximum which is located at Qmax — 0.8 and a location of 
the maximum is shifting from Q ~ 0.7. to Q ~ 0.9 as A increasing. Until reaching Qmax the 
curves are slowly growing and after reaching Qmax they are sharply fiowing down. Fig. 7 shows 
that the field decays faster if it has greater charge. And similar to the spinor case the decay of 
the scalar fields becomes slower as A grows. Both the QN frequencies and the damping rates 
are decreasing with A, although the manner of the decrease is different. 

Fig. 8 displays QNMs for different fixed Q. The curves on Fig. 8 are parametrized by the 
cosmological constant which is changing along the curves. We note that Ima; is related to Rea; 
not linearly in contrast to the spinor field case (Figs. 4, 5). 

Making use of the first order WKB method it is also possible to find the asymptotic formula 
of large / for calculation quasinormal modes lo for the scalar field. 



= Co (/ + 1/2) - z (n + 1/2) + + O (1//) 



(19) 



where 



Cf = TT^o"' {-126QV2 + Q'^rl [378Mro + (-135 + 18Ar2)] + 
o 

rl [-270MVo - 6Mr2 (-30 + 'ikr^) + rjj (-27 + 3Arg)] , 



(20) 
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= -^0 (-eg^ + 2iMgVo - im^rl - 9gV2+ (21) 

IbMrl - 3r^ + 2g2Ar^ - 3MAr^ + Ar^) x 
(42g^ - 12MgVo + + A^Q^l 

-mMrl + 9r^ - eg^Ar^ + 6MAr^ - Arg) . 

In expressions (fTI?j) - (PT|) Cq, tq are the same as those in formulas (fT^ . (fT^ of Sect. 2. By setting 
g = one can make sure that the formula reproduces the corresponding formula of the 
paper [231 obtained for the SdS black hole for the scalar field. We note that the asymptotic 
formula for the RN black hole in asyptotically flat spacetime for the charged scalar field was 
not found till now. The formlula for such a particular case can be obtained by setting A = 
in (Hl-jni). 

4 Conclusions 

We considered the decay of the charged spinor and scalar field near the Reissner-Nordstrom 
black hole in de Sitter spacetime. We calculated the quasinormal frequencies of the fields 
and their damping rates for the lower overtone number. The analysis of the results allows to 
formulate the following conclusions, a) In the spinor case for a fixed cosmological constant A 
the quasinormal frequencies of the charged field are lower and the damping rates are higher 
than those of the neutral field if eQ > and vice versa if eQ < 0; b) In the scalar case for 
a fixed cosmological constant A the fields possessing greater charge have higher quasinormal 
frequencies and decay faster; c) The quasinormal frequencies of both the spinor and scalar fields 
are lower in the spacetimes with greater cosmological constant A. The decay of both the spinor 
and scalar fields is slower in the spacetimes with greater cosmological constant A; d) In the 
case of the charged spinor field the dependence of Imu on Keu is approximately linear on the 
diagram parametrized by A; e) In the case of the charged scalar field the dependence of Imu; 
on Recu is not linear on the diagram parametrized by A; f) The dependence of Imu and Keu 
of the charged spinor field on the black hole charge Q is different for e > and e < for all 
values of A. 

We derived the analytic formula for calculation of u for large values of the multipole number 
both for the spinor and scalar fields. The formula for the spinor field shows that for asymptot- 
ically large multipole number the damping rate does not depend on the charge e of the spinor 
field. 



11 



References 

[1] V. p. Frolov, I.. D. Novikov, Black hole physics: basic concepts and new developments 
(Kluwer Academic publishers, 1998). 

[2] K. Kokkotas and B. Schmidt, Liv. Rev. Rel. 2, 2 (1999); H. R NoUert, Class. Quant. 
Grav.16, R159 (1999). 

[3] S. Hod, Phys. Rev. Lett. 81, 4293 (1998). 

[4] O. Dreyer, Quasinormal Modes, Phys. Rev. Lett. 90, 081301 (2003). 
[5] J. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). 
[6] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998). 

[7] S. Kalyana Rama and B. Sathiapalan, Mod. Phys. Lett. A 14, 2635 (1999). 

G. T. Horowitz and V. E. Hubeny, Phys. Rev. D 62, 024027 (2000). 
[9] V. Cardoso and J. P. S. Lemos, Phys. Rev. D 63, 124015 (2001). 
[10] D. Birmingham, L Sachs, and S. N. Solodukhin, Phys. Rev. Lett. 88, 151301 (2002). 
[11] R. A. Konoplya, Phys. Rev. D 66, 044009 (2002); |hiH^h/0205142 



[12] V. Cardoso, R. Konoplya, and J. P. S. Lemos, Phys. Rev. D 68, 044024 (2003); 
|i?^/0305037| 



[13] R.A. Konoplya, Phys. Rev. D 70, 047503 (2004); |h^h/0406100[ 
[14] S. Perlmutter et al., Astrophys. J. 483, 565 (1997). 

[15] R. R. Caldwell, R. Dave, and P. J. Steindardt, Phys. Rev. Lett. 80, 1582 (1998) 

[16] P. M. Garnavich et al., Astrophys. J. 509, 74 (1998). 

[17] H. Otsuki and T. Futamase, Prog. Theor. Phys. 85, 771 (1991). 

[18] I. G. Moss and J. P. Norman, Class. Quant. Grav. 19, 2323 (2002). 

[19] V. Cardoso and J. P. S. Lemos, Phys. Rev. D 67, 084020 (2003). 

[20] C. Molina, Phys. Rev. D 68, 064007 (2003). 

12 



[21] A. Maassen van den Brink, Phys. Rev. D 68, 047501 (2003). 

[22] R.A. Konoplya, Phys. Rev. D 68, 124017 (2003); |hep-th/0309030l 

[23] A. Zhidenko, Class. Quant. Grav. 21, 273 (2004); gr-qc/03070T2l 

[24] R. A. Konoplya, A. Zhidenko, JHEP 0406, 037 (2004); |hep^h/0402080 

[25] J. Jing, Phys. Rev. D 69, 084009 (2004); |g^/ 0312079 

[26] H. T. Cho, Phys. Rev. D 68, 024003 (2003). 

[27] R. A. Konoplya, Phys. Rev. D 66, 084007 (2002); |g^/0207028 

[28] R.A. Konoplya, Phys. Lett. B 550, 117 (2002); |g^/ 0210105 

[29] H. Zhang and Wei Zhou, Class. Quant. Grav. 21, 917 (2004); |g^/0312029| 

[30] Wei Zhou and Jian-Yang Zhu, Int. J. Mod. Phys. D 13, 1105 (2004); !gr-qc/0309071| 

[31] D. R. Brill and J. A. Wheeler, Rev. Mod. Phys. 29, 465 (1957). 

[32] R. A. Konoplya, Phys. Rev. D 68, 024018 (2003); igr-qc/0303052] 

[33] S. W. Hawking, G. F. R. Ellis, The Large Scale structure of Space Time (Cambridge 
University Press, Cambridge, England, 1973). 

[34] S. Iyer and C. M. Will, Phys. Rev. D 35, 3621 (1987). 



13 



